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Cotunneling transport through quantum dots weakly coupled to non-collinear ly magnetized leads 
is analyzed theoretically by means of the real-time diagrammatic technique. The electric current, dot 
occupations, and dot spin are calculated in the Coulomb blockade regime and for arbitrary magnetic 
configuration of the system. It is shown that an effective exchange field exerted on the dot by 
ferromagnetic leads can significantly modify the transport characteristics in non-collinear magnetic 
configurations, in particular the zero-bias anomaly found recently for antiparallel configuration. For 
asymmetric Anderson model, the exchange field gives rise to precession of the dot spin, which leads 
to a nonmonotonic dependence of the differential conductance and tunnel magnetoresistance on the 
angle between magnetic moments of the leads. An enhanced differential conductance and negative 
TMR are found for certain non-collinear configurations. 
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INTRODUCTION 



Transport through quantum dots coupled to ferro- 
magnetic leads was a subject of extensive considera- 
tions in the last few years. ^^^^^^^ Most of the works con- 
cerned theoretical description of spin-polarized trans- 
port in the strong coupling regime, where the Kondo 
physics emerges, ^i^iii^i^ii^iii as well as in the weak cou- 
pling regime. When the energy of electrons of the source 
lead is larger than the corresponding charging energy of 
the dot, the current in the weak coupling regime flows 
through the system due to sequential (one by one) tun- 
neling via discrete levels of the quantum doti^^ On the 
other hand, if the applied voltage is below threshold for 
sequential tunneling, first-order transport is suppressed 
and the system is in the Coulomb blockade regime. In 
this transport regime the current can still be mediated by 
cotunneling processes, which involve correlated tunneling 
through virtual states of the quantum doti^^^ 

Sequential transport through a single-level quantum 
dot coupled to ferromagnetic leads was studied for 
both collinear— and non-collinear— configu- 
rations of the electrodes' magnetic moments. Spin- 
polarized transport in the cotunneling regime has been 
addressed mainly for collinear systems, ^^^^^^^^ although 
some aspects of cotunneling through quantum dots cou- 
pled to non-collinearly magnetized leads have also been 
considered It has been shown^^ that the interplay 
of single-barrier and double-barrier cotunneling processes 
gives rise to zero-bias anomaly in the differential conduc- 
tance G when the leads' magnetic moments are antipar- 
allel. This, in turn, leads to the corresponding zero-bias 
anomaly in the tunnel magnetoresistance (TMR). More- 
over, the even-odd electron number parity effect in TMR 
has also been found. On the other hand, if the dot 
hosts two orbital levels and the ground state is a triplet, 
zero-bias anomaly develops due to cotunneling through 



singlet and triplet states of the dotM Very recently, ex- 
perimental realizations of quantum dots connected to fer- 
romagnetic leads have also been presented. Such systems 
may be realized for example in granular structures, 
carbon nanotubes,^^^^^^^^'^^ single molecules, magnetic 
tunnel junctions, or single self-assembled InAs quan- 
tum dots.^^ 

The main objective of this paper is to address the prob- 
lem of spin-polarized cotunneling through a quantum 
dot coupled to ferromagnetic leads with non-collinearly 
aligned magnetic moments. The external magnetic field 
required to control magnetic configuration of the sys- 
tem is assumed to be small enough to neglect the Zee- 
man splitting of the dot level. When the system is 
in a non-collinear magnetic configuration, an effective 
exchange field can significantly modify the transport 
characteristics We show that in the case of symmet- 
ric Anderson model, e = where £ is the dot level 
energy and U denotes the Coulomb correlation parame- 
ter, the zero-bias anomaly in differential conductance de- 
creases monotonically with increasing deviation of mag- 
netic configuration from the antiparallel alignment, and 
eventually vanishes in the parallel configuration. On the 
other hand, the effect of exchange field is most visible 
when the system is described by an asymmetric Ander- 
son model, e 7^ leading to a nontrivial variation 
of the differential conductance with the angle between 
magnetizations of the leads. 

The anomalous behavior of the differential conduc- 
tance leads to a related anomaly in TMR,^^ which for 
arbitrarily aligned magnetizations of the leads can be 
defined as TMR = [/p - I{lp)]/ I{lp), where /p is the 
current flowing through the system in the parallel con- 
figuration {(f> = 0) and I{(p) is the current flowing in the 
non-collinear magnetic configuration (i.e. when there is 
an angle (p between the leads' magnetic moments). We 
show that the angular dependence of TMR strongly de- 
pends on the system parameters - in the case of asym- 
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metric Anderson model TMR can become negative when 
the leads are aligned non-collinearly. 

The paper is organized in the following way. In sec- 
tion [ni we present the model of single-level quantum 
dot coupled to non-collinearly magnetized leads. In sec- 
tion mil we briefly describe the real-time diagrammatic 
technique used to calculate transport properties. We 
also present the perturbation expansion in the Coulomb 
blockade regime for the case of arbitrary magnetic con- 
figuration of the system. Numerical results followed by 
discussion are presented in section HVl while the conclu- 
sions are given in section IVl 
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II. MODEL 

The schematic of the system considered is shown in 
Fig. [TJ The system consists of a single- level quantum dot 
coupled to ferromagnetic leads, whose spin moments, Sl 
and Sr for left and right lead, form an arbitrary con- 
figuration. The magnetic configuration of the system is 
described by the angle - when 99 = = tt) the 
system is in the parallel (antiparallel) magnetic configu- 
ration, see Fig. [TJ 

The system is described by an Anderson-like Hamilto- 
nian of the form 



(1) 



where the first two components describe noninteracting 
itinerant electrons in the leads. Ha = Xl/ccr ^c^kacl^ka^o^ka 
for the left {a = L) and right {a = R) lead in the lo- 
cal reference frames. The energy of an electron with 
spin a and wave number k in the lead a is denoted 
by Saka^ while cj^^^ {caka) IS the corresponding creation 
(annihilation) operator. The dot is described by i^D, 
which in the global reference frame can be expressed 
as Hd — Xl^^i I eS^dcr + Un^ni^ where e is the single- 
particle energy, while dj. {da-) creates (destroys) a spin-a 
electron in the dot. The second term of the dot Hamilto- 
nian describes the Coulomb interaction of two electrons 
of opposite spins residing in the dot, with U denoting the 
corresponding correlation energy. For the system geom- 
etry shown in Fig. [H the tunnel Hamiltonian is given by 
i^T = ^TL + ^TR , with 
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(2) 



COS - -TR_c{^fc_sin-jdt 



while -ffxR reads 

k 

+ (Tr+4^^ sin I + T^-al^k- cos |) + h.c] (3) 

where T^a denotes the tunnel matrix elements between 
the dot states and the majority (a = +) or minority 
(cr = — ) electron states in the lead a when = 0. 



FIG. 1: (color online) Schematic of a quantum dot coupled to 
ferromagnetic leads with non-collinearly aligned magnetiza- 
tions. The net spin moments of the left Sl and right Sr lead 
form an angle f. There is a symmetric bias voltage applied 
to the system. 



The coupling of the dot to ferromagnetic leads in the 
case of collinear magnetic configuration is described by 
F^ = 27:\Tacr\^ pacri whcrc p^cr is the density of states 
for the majority (a = +) and minority (a = — ) elec- 
trons in the lead a. It is convenient to express the cou- 
pling parameters in terms of spin polarization defined 
as Pcx = (r+ - F^)/(r+ -h F^). Thus, the coupling 
strength can be written as F^ = Fq,(1 ± Pq,), where 
Fq, = (F+ + F~)/2. In the following, we assume that the 
couplings are symmetric, Fl = Fr = F/2. As reported 
in Ref. 38, in the weak coupling regime typical values of 
the dot-lead coupling strength F are of the order of tens 
of /ieV. 

We assume that the magnetic configuration of the sys- 
tem can be changed by applying a weak external mag- 
netic field. The corresponding Zeeman energy is then 
very small and can be neglected. Thus, following pre- 
vious theoretical calculations,^^ we assume that the dot 
energy level is spin-degenerate. In the case of quantum 
dots, the magnetic fields that produce experimentally 
measurable effects are rather high, ~ 10^ OeM 



III. METHOD AND TRANSPORT EQUATIONS 

In order to calculate transport properties of the system 
we employ the real-time diagrammatic technique. ^^^^^^ 
This technique is based on a systematic perturbation ex- 
pansion of the dot density matrix in the dot-lead coupling 
strength, which allows one to calculate transport charac- 
teristics order by order in tunneling processes. In the 
following, we address the problem of spin-polarized tun- 
neling in the Coulomb blockade regime. In this trans- 
port regime the first-order (sequential) tunneling pro- 
cesses are exponentially suppressed due to the charging 
energy and the current is mediated by second- (cotunnel- 
ing) or higher-order tunneling processes.— li^ii^ 

We note that although the approach based on the 
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second-order perturbation theory and standard rate 
equation would lead to reasonable results for the 
Coulomb blockade regime and collinear magnetic config- 
urations, it would be insufficient in the case of arbitrarily 
aligned magnetizations of the leads. This is because this 
approach does not allow for the effects of exchange field. 
On the other hand, the advantage of the real-time di- 
agrammatic technique is that it takes into account the 
exchange field in a fully systematic way. 

In our considerations we assume that the correlation 
parameter U between two electrons in the dot is finite, 
which results in four dot states |x), with x = O^T^i^d, 
respectively for the empty dot, singly occupied dot with 
spin-up or spin-down electron, and for two electrons in 
the dot. The density matrix of the dot involves then six 
elements. 



whereas the second-order master equation reads 



/PO \ 

P/ P/ 

P/ P/ 

V OPiJ 



(4) 



The diagonal elements of the density matrix correspond 
to the respective occupation probabilities, while the off- 
diagonal elements Pj and Pj^ describe the dot spin 



with S:, = ReP/, Sy = ImP/, and S, = (p/ - P/) /2. 

The time evolution of the dot density matrix is described 
by the Dyson equation, which can be transformed into a 
master-like equation. The master equation allows one 
to determine the elements of the density matrix. In the 
steady state and for spin-degenerate dot level it can be 



written a; 



,39,40 
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From the above equations and using the normalization, 
^x^x^^^ = Sj^o^ one can find the density matrix ele- 
ments in the zeroth and first order. Thus, the key point 
is to calculate the first- and second-order self-energies. 
This can be done with the aid of the following rules in 
the energy space: 

1. Draw all topologically different diagrams with fixed 
time ordering and position of vertices. Connect the 
vertices by tunneling lines. Assign the energies of 
respective quantum dot states to the propagators 
and a frequency uo to each tunneling line. 

2. Tunneling lines acquire arrows indicating whether 
an electron leaves or enters the dot. For tunneling 
lines changing the spin from a to a' assign 7^^^ (cj) 
ha^^ (cc;)] if the tunneling line goes backward (for- 
ward) with respect to the Keldysh contour and elec- 
tron tunnels between the lead a and the dot. 

3. For each time interval on the real axis limited by 
two adjacent vertices assign a resolvent 1/{AE + 
iO+), with AE being the difference of all energies 
going from right minus all energies going from left. 

4. Each diagram gets a prefactor (—1)^, with 7 being 
the number of internal vertices lying on the back- 
ward propagator plus the number of crossings of 
the tunneling lines, plus the number of vertices con- 
necting state |d) with state | |). 

5. Integrate over all frequencies and sum up over the 
reservoirs. 



where E^,^^^ is the irreducible self-energy corresponding 

X2X2 

to transition forward in time from state |xi) to |xi) and 
then backward in time from state |x2) to \X2)- 



A. Dot occupations and spin 

We consider transport in the Coulomb blockade 
regime, |£ + /7| F, /cbP, and when the dot is singly 
occupied, e < < £ -\- U . In order to calculate the 
dot occupations and spin, we expand the self-energies. 
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E. T Ti^j^^ , and density matrix elements, 
J=i X2X2 ' ^ ' 

^^•^Q P^^ , order by order in tunneling pro- 
cesses. As a consequence, the master equation in the 
first order is given by 



0= y p^/^'^E^^ 

Z-^ X2 X2 



XiXi 

X2 



(1) 



(6) 



The factor 7^^^'(cj) is defined as, it'^'^'i^) 
r--7±(c^-/i«), with 



— sin(^(F+-Fj^) 



^ R 
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F+ cos^ f + ^i^^ f 

- • 2 ^ , -P- 2 ^ 
^ . ^ ^. sm — + I p cos — 
27r V ^ 2 ^ 2 



(8) 
(9) 
(10) 



while Fl^ is given by changing R ^ L and setting ip = 
0. Here, is the Fermi-Dirac distribution function, 
and /~ = 1 — /+. We note that for arbitrary SU(2) 
rotations F^^ = (F^^) , while for rotations around the y 
axis (considered in this paper) F^^ = TjJ^. 

By using the above rules, we find some useful general 
relations for the self-energies which hold in the first and 
second order 
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with X = O5 T, d and j = 1, 2. In addition, 



(i) 



(12) 



(13) 



for x,x' = cr,cr (x 7^ xO^ and S^^^ ^ = 0, as transition 
between such states has to involve at least two tunneling 
lines. In general, the self-energies are complex. However, 
as we show in the following, in the case of the Coulomb 
blockade regime one only needs to determine the imag- 
inary parts of the second-order self-energies. This con- 
siderably simplifies numerical calculations, since the self- 
energies can be determined analytically. 

First of all, we note that in the Coulomb blockade 
regime the first-order self-energies associated with en- 
ergetically forbidden transitions are exponentially small, 
i.e., E-^(i) = 0, for x = 0, T, I, d, and S^^^'^ = 0, for x = 
0,d. Furthermore, one also has ImS^^^"*^^ = ImS^^^"^^ = 
ImE^^^^^ =0. As a consequence, from the first-order 



master equation, Eq. (|6]), one finds, Pq 



(0) 



and 5*^^^ = 0, whereas the occupations P^' ' P, 
are related by the equation 



t(0) ^|(o) 



= 0, 
and 



P/'°^s|f^^P/^°^E|f^^5(°)Elf^^=0. (14) 

This equation, however, cannot be solved, as it involves 

three unknown variables. To solve for , , and 

Si'\ it is necessary to include the second-order pro- 
cesses. If the system is in a collinear configuration (either 
parallel or antiparallel) the first two self-energies vanish 

m^^^ = Sjj^^^ = 0, and Eq. ^ yields Si^^ = 0. It is 
also worth noting that the above equation describes the 



effect of an exchange field which is exerted by ferromag- 
netic leads with non-collinearly aligned magnetizations.^^- 
This exchange field results mainly from the first-order 
processes. To realize the meaning of Eq. (p!4|) we rewrite 
this equation in a more explicit form (see the formulas 
given in Appendix A) 



4^ 



, (15) 



where 



Br =rrRe 



27rkBT J 

U - lla 



(16) 



describe the exchange field effects, with ^ denoting the 
digamma function. The exchange field results from the 
spin-dependent virtual processes between the dot and 
leads. The angular variation of the exchange field is 



mainly governed by B^i ~ pRrRsin(/:? {B{; = in the 
assumed coordinate system). We note that B^^ vanishes 
for collinear magnetic configurations as well as for non- 
magnetic leads. As follows from Eq. ([15]), the exchange 
field in the first order gives rise to the rotation of the dot 
spin in the x — z plane^i Furthermore, we note that it is 
a purely many-body effect, as the exchange field vanishes 
in the noninteracting case, t/ = 0. 

Taking into account the aforementioned results and re- 
marks, one can set up the second-order master equation, 
Eq. (O, for the Coulomb blockade regime. From the 
master equation in the second-order one gets a set of 
equations which, together with Eq. ([El), allows one to 



P^^^^ and P 



d(i) 



and also 



S^^ and The equations which enable full solution 

of the problem can be written in a matrix form as 
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(17) 



where the last row is due to the normalization condi- 

|(0) I (0) 

tion, P| + ^1 =1- Equation (|17j) involves only the 
imaginary parts of the second-order self-energies which 
can be calculated analytically. The formulas for the self- 
energies appearing in Eq. ([T7|) and having at the ends the 



off-diagonal states are given in Appendix A. 

From the second-order master equation one also ob- 

t(i) 1(1) 

tains an additional equation which involves PJ , P| , 
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+Pj^'^E\f^ +Pl^'^E\f^ +Si'^E]f^ . (18) 



This equation is analogous to Eq. ([T4j) in the first order. 

To solve it for , P/ , and si^^ one would need 

to include the higher-order (third-order) diagrams. In 
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the case of collinear magnetic configurations, S 
Sjj^^^ = for j = 1,2, and from Eq. ([18]) one finds 

Si^'' = 0, while S^-' = from the first order master 
equation. Thus, if the magnetizations of the leads are 
aligned collinearly, Eq. (pT|) reduces to the one presented 
in Ref. M- 



B. Green's functions and the current 

The current flowing through the system can be ex- 
pressed in terms of the Green's functions as^ 



with 



^ I du;Tr{[rL/+H-rR/+H] 



[G> (c.) - G< (c.)] + (Fl - Tr) G< {uj) } (19) 



case of Coulomb blockade regime, the first-order current 
is exponentially suppressed and the dominant contribu- 
tion comes from the second-order processes. Therefore, 
one only needs to calculate the Green's functions in the 
second order. This can be done using the following rules 
in the energy space: 

1. Assign the corresponding external and, in higher 
orders, internal vertices to a bare Keldysh contour. 
Connect the external vertices by a virtual line and 
the internal vertices by tunneling lines. Assign the 
energies of the respective quantum dot states to 
the propagators and a frequency uj to each tunnel- 
ing line and the virtual line. Multiply the diagram 
with an element of the reduced density matrix cor- 
responding to the initial state. 

2. See rule (2) in calculation of self-energies. 

3. See rule (3) in calculation of self-energies. 

4. See rule (4) in calculation of self-energies. 

5. Integrate over all frequencies except for the virtual 
line and sum up over the reservoirs. 

For the lesser Green's function an additional minus sign 
appears, which is due to reversed direction of the vir- 
tual line. The formulas for the Green's functions in the 
Coulomb blockade regime are given in Appendix B. 



(20) 



rR = rR( ^ ' ^^---^ ^^--^ ) , (21) 

PR sm ip 1 — PR cos (fi I 



The coupling matrices are given by 

1 + PR cos PR sin if 



for the right lead and 



rL = rL 



1+Pl 

I-PL 



(22) 



for the left lead. The Green's functions G^^,{lj) are de- 
fined as Fourier transforms of 



G>^,it) = -i{d,it)di,m ■ 



(23) 
(24) 



In order to calculate the current in each order we ex- 
pand the Green's functions systematically in the coupling 
strength F— 

oo oo m 

Gt, = Y: GflT^ = ^ ^ Gfi?'™-") . (25) 



m=0 



m=0 n=0 



Thus, for the first order one has Cji?^ = whereas 

' (T(T (7(7 ' 

in the second order with respect to F the Green's func- 
tions are given by Gf[]^ = G^i?'^^ + Gf^}'^\ In the 

c> J (7(7' (7(7' (7(7' 



IV. RESULTS AND DISCUSSION 

In the following we discuss numerical results on elec- 
tronic transport through a single-level quantum dot cou- 
pled to non-collinearly magnetized leads. We first present 
the results for a symmetric Anderson model, e = 
then we proceed with the discussion of an asymmetric 
Anderson model, e ^ —U/2. If e = — 1//2, the ex- 
change field vanishes and both differential conductance 
and TMR change monotonically when going from paral- 
lel to antiparallel magnetic configurations. However, for 
e 7^ the effects of exchange field become important 

and lead to nontrivial behavior of transport characteris- 
tics. The symmetric and asymmetric Anderson models 
can be realized experimentally by applying a gate voltage 
to the dot. With the gate voltage, one can continuously 
shift position of the dot level (while U is unaffected) and 
thus study the cross-over from symmetric to asymmetric 
models. 



A. Symmetric Anderson model 

The differential conductance G as a function of the 
bias voltage and the angle is shown in Fig. [2ja), while 
Fig. El^b) displays the bias dependence of G for several 
values of ip. In the collinear, i.e. parallel and antiparallel 
magnetic configurations we recover the earlier results 
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FIG. 2: (color online) (a) Differential conductance G ■ 
dl/dV as a function of the bias voltage V and the ang^ 
(f for the symmetric Anderson model. The parameters ar( 
/cbT O.Sr, s -15r, [/ 30r, and pL = Pr = p = 0.1 
Part (b) shows the bias dependence of differential conductance 
for several angles. 



with the zero-bias anomaly present in the antiparallel 
configuration and a parabolic dependence of G on bias 
voltage for the parallel configuration. Now, we find a 
monotonic variation of the anomaly with the angle be- 
tween leads' magnetic moments. When the angle in- 
creases from zero to tt, the anomaly emerges at small 
values of (p and its relative height increases with increas- 
ing angle, reaching a maximum value at (f = tt. The 
anomaly results from the interplay of nonequilibrium spin 
accumulation in the dot and the competition between the 
single-barrier and double-barrier cotunneling processes. 
The physical mechanism of this anomalous behavior was 
discussed in detail in Ref. [23], therefore, we will not de- 
scribe it here. Instead, we will focus on its variation with 
the angle between magnetic moments of the leads and on 
the difference between symmetric and asymmetric An- 
derson models. 

In the linear response regime and for symmetric An- 
derson model the exchange field, Eq. ([IS]), is negligible 
and the average dot spin tends to zero. It is therefore 
relatively easy to find an approximate formula for the 
differential conductance. For e = Pl = Pr = 

and at low temperatures, the angular dependence of the 
linear conductance is given by 



2p2 r 



G = 



3-/ (l + 2sin2|) 



(26) 
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FIG. 3: (color online) (a) Tunnel magnetoresistance as a func- 
tion of the bias voltage V and the angle cp for the symmetric 
Anderson model. The parameters the same as in Fig. [2] Part 
(b) shows the TMR as a function of the bias voltage for several 
values of cp. 



The variation of G with (f at low bias is thus character- 
ized by the factor 1 + 2 sin^((/:)/2), which leads to maxi- 
mum (minimum) conductance in the parallel (antiparal- 
lel) magnetic configuration. Such behavior is typical of a 
normal spin- valve effect. 

The TMR ratio, defined in the introduction, is shown 
in Fig. O where part (a) displays the angular and bias 
dependence, while part (b) shows the bias dependence 
of TMR for several values of the angle cp. The zero- 
bias anomaly in the differential conductance (see Fig. [2]) 
leads to the corresponding anomaly (dip) in the TMR 
at small bias voltages. For collinear configurations we 
recover again the results presented in Ref. [2J]. The dip 
in TMR decreases when magnetic configuration departs 
from the antiparallel one, and eventually disappears in 
the parallel configuration. The variation of TMR with 
the angle is monotonic, similarly as the angular variation 
of the differential conductance. 

With the same assumptions as those made when de- 
riving Eq. ([26]) , one can approximate the dependence of 
TMR on the angle cp at zero bias by 



TMR: 



2p' 



^sin^f 



3-p2(i + 2sin2f) 



(27) 



Now, the angular dependence of TMR is governed by 
sin^((/?/2), which gives maximum TMR in the antiparallel 
configuration and zero TMR in the parallel one. 
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FIG. 4: (color online) The Sx (a) and Sz (b) components 
of the dot spin as a function of the bias voltage for several 
values of the angle cp as indicated in the figure. The other 
parameters are the same as in Fig. (2] In part (a) the curves 
for (f = and cp = tt coincide. 



tributions to the exchange field do not cancel then, and 
the resulting effective exchange field becomes nonzero. 
This has a significant impact on transport properties, as 
we show and discuss below. 

First, we note that the leading contribution to the ex- 
change field comes from the first-order diagrams, whereas 
the current flows due to the second-order tunneling pro- 
cesses. This leads to two different time scales which de- 
termine transport characteristics, as will be discussed be- 
low in more details. Transport properties are thus a re- 
sult of the interplay between the first- and second-order 
processes. We emphasize that although the first-order 
processes do not contribute directly to the current, they 
influence transport via modification of the dot spin. 

The strength of effective exchange field is determined 
by deviation of the Anderson model from the symmetric 
one, described quantitatively hy 2e-\-U (with 2e-\-U = 
for the symmetric model). On the other hand, the cotun- 
neling processes contributing to the current are propor- 
tional to the bias voltage. Thus, to parameterize relative 
magnitude of the two processes we define a dimensionless 
parameter x = {2e + U)/\eV\. 

In Fig. [5^ we show the differential conductance as a 
function of angle and bias voltage for an asymmetric An- 
derson model, £ 7^ — [//2, while Fig. [5)3 displays the bias 
dependence of G for different values of cp. As one can 

see ^"h<=» loA^f-KiftQ rlifF<=»r<=»nf ial nnnrlnnf anp^i K<=»r'nm<=»« f^Jl- 



In Fig. [4] we show the average components Sx and Sz 
of the dot spin as a function of the bias voltage, calcu- 
lated for several values of the angle cp. As one can see, 
Sx = and = in the parallel magnetic configu- 
ration {(f = 0), while in the antiparallel configuration 
{(p = tt) Sx = whereas Sz is finite. The y compo- 
nent of the dot spin vanishes, Sy = 0^ irrespective of 
magnetic configuration of the system, so it is not shown 
in Fig. m As a consequence, up to the second order in 
tunneling processes, the average dot spin is in the x — z 
plane, which is a characteristic feature of the symmetric 
Anderson model, e = —U/2. In other words, there is 
no precession of the spin out of the x — z plane, that is 
out of the plane formed by magnetic moments of the two 
electrodes. This is due to the cancellation of different 
contributions to the exchange field in the case of sym- 
metric Anderson model. The absence of exchange field is 
responsible for the monotonic angular variation of both 
the differential conductance and TMR, discussed above, 
and also for vanishing of Sy component of the dot spin. 



B. Asymmetric Anderson model 

The transport characteristics described above for the 
symmetric Anderson model are significantly modified 
when the model becomes asymmetric, i.e. when e ^ 
—U/2. This can be realized by shifting the dot level po- 
sition by a gate voltage applied to the dot. Different con- 




FIG. 5: (color online) (a) The differential conductance G — 
dl/dV as a function of the bias voltage V and the angle cp for 
the asymmetric Anderson model e — — 12r and U — SOP. The 
other parameters are the same as in Fig. [2] Part (b) shows the 
bias dependence of differential conductance for several values 
of the angle (p. 
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FIG. 6: (color online) (a) Tunnel magnetoresistance as a func- 
tion of the bias voltage V and the angle (p for the asymmetric 
Anderson model. The parameters the same as in Fig. [5] Part 
(b) shows the TMR as a function of the bias voltage for several 
values of cp. 



FIG. 7: (color online) The differential conductance G (a) and 
tunnel magnetoresistance (b) as a function of the angle cp for 
several values of x = {2e + U)/\eV\ and for e = — 12r. The 
other parameters are the same as in Fig. (2] 



hanced in a certain range of the angle (p^ leading to a 
nonmonotonic dependence of the conductance G on the 
angle between the leads' magnetizations. This nonmono- 
tonic behavior is also visible in TMR. The density plot 
of TMR is shown in Fig. [6^1, whereas the TMR as a func- 
tion of bias voltage for different values of (p is shown in 
Fig. [6)3. There is a range of the angle cp between mag- 
netic moments of the leads, where TMR changes sign and 
becomes negative in the small bias region, i.e. the corre- 
sponding conductance is larger than that in the parallel 
configuration. To discuss and account for the influence 
of exchange field, we further present one-dimensional fig- 
ures displaying explicitly the angular dependence of both 
G and TMR. 

The angular variation of the differential conductance 
and TMR for an asymmetric Anderson model is shown 
in Fig. [71 for several values of the parameter x, x = 
{2£ + U)/\eV\. In this figure both e and U are kept con- 
stant {e = —2U /5 and U = SOF), so the variation of x is 
only due to the change of the bias voltage. When increas- 
ing X (decreasing bias voltage), one increases the relative 
role of the first-order processes (leading to the spin pre- 
cession), which considerably modifies the angular depen- 
dence of the differential conductance and TMR. As one 
can see in Fig. [7^1, the differential conductance for x 1 
does not change monotonically with the angle cp. It in- 
creases first with increasing (p^ reaches a local maximum, 
and then drops to a local minimum for (p close to the an- 



tiparallel configuration, and again slightly increases when 
approaching (p = tt. Thus, the maximum and minimum 
differential conductance occur for non-collinear configu- 
rations, and not for parallel and antiparallel ones. The 
local maximum in G at the antiparallel configuration be- 
comes smaller when the effect of exchange field is sup- 
pressed (x decreases) and eventually disappears for x = 1 
by transforming into the global minimum of G in the an- 
tiparallel configuration. Similarly, the maximum at a 
non-collinear configuration disappears with decreasing x 
and changes into a global maximum in the parallel con- 
figuration. 

The most characteristic features of transport charac- 
teristics in the presence of exchange field are the en- 
hanced differential conductance at a non-collinear align- 
ment and its rapid drop when approaching the antipar- 
allel configuration. This is most visible in the curve for 
X = 50 in Fig. [7^1. In this case the differential conduc- 
tance increases until the configuration becomes close to 
the antiparallel one, and then drops rapidly to the value 
of G in the antiparallel configuration. The key role in 
this behavior is played by the first-order processes giving 
rise to the exchange field. These processes lead to the 
precession of spin in the dot, which facilitates tunneling 
processes and leads to an increase in the conductance as 
compared to the parallel configuration. When the con- 
figuration becomes close to the antiparallel one, the first- 
order processes become suppressed and the conductance 
drops to that for antiparallel alignment. For smaller x the 
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FIG. 8: (color online) The differential conductance (a) and 
tunnel magnetoresistance (b) as a function of angle between 
the leads' magnetizations for several values of the level posi- 
tion £, as indicated in the figure, and for 1/ = F. The other 
parameters are the same as in Fig. [2] 



drop in conductance begins at a smaller angle, which fol- 
lows from the fact that now the relative role of exchange 
field decreases due to increased cotunneling rates. 

The nonmonotonic behavior of the differential conduc- 
tance with (f leads to a nonmonotonic dependence of 
TMR. This is shown in Fig. [7)3 for several values of x. 
The effects due to exchange field give rise to a local min- 
imum in TMR at a non-collinear magnetic configuration. 
Moreover, in this transport regime TMR changes sign 
and becomes negative. When the magnetic configuration 
is close to the antiparallel one, TMR starts to increase 
rapidly reaching maximum for (p = it. The negative TMR 
and its sudden increase when the configuration tends to 
the antiparallel one are a consequence of the processes 
leading to nonmonotonic behavior of the differential con- 
ductance, as described above. 

To understand more intuitively the above presented 
behavior of the differential conductance and TMR at low 
bias voltage and close to the antiparallel configuration, 
one should consider two different time scales. One time 
scale, Tprec, IS established by the virtual first-order pro- 
cesses responsible for the spin precession due to exchange 
field. 



2h 



p sin cp In 



(28) 



The second time scale, Tcot, is associated with second- 
order processes which drive the current through the sys- 
tem. At low temperature and bias voltage, the cotunnel- 



ing rate can be expressed as 
F2 



^cot 



^^(l+p)(l-,cos,)^l^. (29) 



The rate Tcot depends linearly on the applied voltage, 
whereas Tprec is rather independent of V. As a conse- 
quence, at low bias and for non-collinear configuration, 
the exchange field plays an important role leading to a 
nonmonotonic dependence of differential conductance on 
the angle between the leads' magnetizations. When mag- 
netic configuration is close to the antiparallel one, the 
spin precession rate is deceased (|Tp^ec| smcp) and, at 
certain angle, the rate of spin precession becomes compa- 
rable to the cotunneling rate. This gives rise to a sudden 
drop (increase) in differential conductance (TMR), as can 
be seen in Fig. [7^(b). We note that the nonmonotonic 
dependence of differential conductance and magnetore- 
sistance has also been observed in quantum dots in the 
strong coupled limit. ^ 

As follows from Eqs. ([28|) and ([29|) . in the low bias volt- 
age regime the influence of the exchange field is mainly 
determined by the magnitude of departure from the sym- 
metric Anderson model, |Tp^ec| ^ In [^/(e + [/)|. The an- 
gular variation of the differential conductance and TMR 
for several values of the level position is shown in Fig. [8l 
When raising the position of the dot level one goes from 
the symmetric to asymmetric Anderson model, increas- 
ing the effect of exchange field, which is most visible for 
e = — lOF. The influence of exchange field is the same as 
in the case of Fig. [71 discussed above. Figure [8] demon- 
strates that at a certain non-collinear magnetic configu- 
ration, by sweeping the gate voltage one can change the 
differential conductance and also tune the TMR between 
the positive and negative values. This can be of some 
importance for future applications in spintronics. 

In the case of symmetric Anderson model the dot spin 
is located in the plane defined by the magnetizations of 
the leads, i.e. in the x — z plane. On the other hand, 
in the case of asymmetric Anderson model, the exchange 
field gives rise to the spin precession in the dot, leading 
to a nonzero Sy component of the dot spin. The com- 
ponents of the spin in the dot are shown in Fig. [9] for 
different position of the dot level. If we start from a 
symmetric Anderson model, then the conductance and 
TMR are independent of the sign of the dot level shift 
(shifting the dot level up or down by the same amount 
leads to equal conductance and TMR). The spin in the 
dot, however, depends on the sign of the shift. This is 
due to the fact that changing sign of the shift is associ- 
ated with a sign change of the exchange field, which is 
proportional to ln\e/ {e -\- U)\. As can be seen in Fig. [9l 
the sign of the exchange field determines the sign of the 
spin components. It is also worth noting that the two 
components, Sx and Sy, disappear in both collinear con- 
figurations (parallel and antiparallel). In turn the Sz 
component vanishes only in the parallel configuration, 
while it remains nonzero in the antiparallel configuration 
due to spin accumulation^ 
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FIG. 9: (color online) The Sx (a), Sy (b), and Sz (c) com- 
ponents of the dot spin as a function of angle cp for several 
values of e and for V = F. The other parameters are the same 
as in Fig. [2] 



defined by the magnetic moments of the leads. 

When the system is described by an asymmetric An- 
derson model, the angular dependence of the differential 
conductance and TMR becomes more complex. First, 
the maximum and minimum values of differential conduc- 
tance occur not for collinear configurations, but when the 
leads' magnetizations are non-collinearly aligned. Sec- 
ond, the conductance is enhanced in a broad range of 
the angle between the magnetic moments of the leads, 
and drops to the minimum value when the configura- 
tion becomes close to the antiparallel one. The conduc- 
tance drop is particularly pronounced for small bias volt- 
ages. This behavior results from the interplay between 
the first-order processes giving rise to the exchange field 
and the second-order processes that drive the current. 
The exchange field also leads to a nonmonotonic angular 
dependence of TMR. For magnetic configurations where 
differential conductance is enhanced as compared to that 
in parallel configuration, the TMR becomes negative. 
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APPENDIX A: SELF-ENERGIES IN THE 
COULOMB BLOCKADE REGIME 



V. SUMMARY AND CONCLUDING REMARKS 

We have considered the cotunneling transport through 
a single-level quantum dot coupled to ferromagnetic leads 
with non-collinearly aligned magnetic moments. Trans- 
port properties have been calculated with the aid of 
the real-time diagrammatic technique which has been 
adapted to the case of deep Coulomb blockade regime. 
The advantage of the real-time diagrammatic technique 
is that it takes into account the exchange field resulting 
from virtual processes between the dot and ferromagnetic 
leads in a fully systematic way. 

In the case of symmetric Anderson model, we have 
found a monotonic variation of the differential conduc- 
tance and TMR with the angle between magnetic mo- 
ments of the leads. This typical angular dependence is 
due to the fact that different contributions to the ex- 
change field cancel each other and there is no spin pre- 
cession of the dot spin - the dot spin remains in the plane 



In this appendix we present the analytical formu- 
las for the first- and second-order self-energies involv- 
ing at the ends the off-diagonal states that are neces- 
sary to determine the dot occupations and spin from 
Eq. (pT|) . Our results are valid in the case of deep 
Coulomb blockade, \s\,\£ + /7| > kBT.T. For the 

first order we find, E^^^'^ = ^fj^^ = 27ri^^r-f, 
^^^(1) ^ ^^^(1) ^ ^JB?^-C^^), and = 
Ea (B^a - B'f^ + Cff - Cf,^), where a = L, R. 

The imaginary parts of the second-order self-energies 
involving off-diagonal states are given by 



(2) 



Q!,/3 a' 



2rf {B, 



2/3 



r) -r)<J<J.<J<7 



2c: 
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the case of the Coulomb blockade regime. The Green's 

functions G^^l'^^ are given by Gai"*^'^^ = —27riPQ^^'^S{u; — 

e), G^i''^^ = 27TiPf^^S{co^- e - U), while Gfi''^^ = 0. 
On the other hand, for G^i?'^^ we find 

■ (7(7 



^>(0,1) 



-27rz 



^1 ^ 



(0) 



+25f [62-r + B2l<5(c.-£) 
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(&rr-^=rr)|, (B2) 



ImS?°'')=27r^(r-+rf)B0, 

a,/3 



(A5) (^<(0,1) 



where a, = L, R and a', <j'' = a, a. The corresponding 
parameters are defined as 



■^na ~ ^a ^n(^ Ma) •> 



(A6) 



^nal3 



= r 



(7(7 T^cr cr 



[^n(^: - /ic) - Xn{e - /i/3)] , (A7) 



with 

^n+i{x) = ^Re 



27rkBT 



In 



27rA:BT 
(A8) 

where is the cutoff parameter and /b denotes the Bose 
function. The parameters G^^ and C^^q ^ are given 



by Eqs. (|A6p and (|A7p . respectively, with e replaced by 
e + U, while AZf^'" = B^f-'" - CZf-'" . The 
other self-energies can be found using the rules and rela- 
tions given in section UlI Ai 



APPENDIX B: GREEN'S FUNCTIONS IN THE 
COULOMB BLOCKADE REGIME 

In the following we present the explicit formulas for 
the Green's functions needed to calculate the current in 



27rz 
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u 



pcr(O) I p, 



with 



L±(7a' 



(CO) 



(B5) 



and c^^^' is given by Eq. ()B5p with e ^ s -\- while 
the other parameters are defined in Appendix A. In the 
above formulas we have cancelled all the terms whose 
contribution to the current is exponentially suppressed, 
i.e., the terms multiplied with S{u; — e) and S{u; — e — U) 

in G^^?'"*^^ and G^^,'^\ respectively, as well as the terms 
multiplied with 5'{(jO — e) and d'ioo — e — U)^ where 5 is 
the Dirac function and 5' its derivative. 
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